New insights into Trypanosoma cruzi genetic diversity, and its influence on parasite biology and clinical outcomes.

Chagas disease, caused by Trypanosoma cruzi, remains a serious public health problem worldwide. The parasite was subdivided into six distinct genetic groups, called "discrete typing units" (DTUs), from TcI to TcVI. Several studies have indicated that the heterogeneity of T. cruzi species directly affects the diversity of clinical manifestations of Chagas disease, control, diagnosis performance, and susceptibility to treatment. Thus, this review aims to describe how T. cruzi genetic diversity influences the biology of the parasite and/or clinical parameters in humans. Regarding the geographic dispersion of T. cruzi, evident differences were observed in the distribution of DTUs in distinct areas. For example, TcII is the main DTU detected in Brazilian patients from the central and southeastern regions, where there are also registers of TcVI as a secondary T. cruzi DTU. An important aspect observed in previous studies is that the genetic variability of T. cruzi can impact parasite infectivity, reproduction, and differentiation in the vectors. It has been proposed that T. cruzi DTU influences the host immune response and affects disease progression. Genetic aspects of the parasite play an important role in determining which host tissues will be infected, thus heavily influencing Chagas disease's pathogenesis. Several teams have investigated the correlation between T. cruzi DTU and the reactivation of Chagas disease. In agreement with these data, it is reasonable to suppose that the immunological condition of the patient, whether or not associated with the reactivation of the T. cruzi infection and the parasite strain, may have an important role in the pathogenesis of Chagas disease. In this context, understanding the genetics of T. cruzi and its biological and clinical implications will provide new knowledge that may contribute to additional strategies in the diagnosis and clinical outcome follow-up of patients with Chagas disease, in addition to the reactivation of immunocompromised patients infected with T. cruzi.

Manipulating multi-level selection in a fungal entomopathogen reveals social conflicts and a method for improving biocontrol traits.

Changes in parasite virulence are commonly expected to lead to trade-offs in other life history traits that can affect fitness. Understanding these trade-offs is particularly important if we want to manipulate the virulence of microbial biological control agents. Theoretically, selection across different spatial scales, i.e. between- and within-hosts, shapes these trade-offs. However, trade-offs are also dependent on parasite biology. Despite their applied importance the evolution of virulence in fungal parasites is poorly understood: virulence can be unstable in culture and commonly fails to increase in simple passage experiments. We hypothesized that manipulating selection intensity at different scales would reveal virulence trade-offs in a fungal pathogen of aphids, Akanthomyces muscarius. Starting with a genetically diverse stock we selected for speed of kill, parasite yield or infectivity by manipulating competition within and between hosts and between-populations of hosts over 7 rounds of infection. We characterized ancestral and evolved lineages by whole genome sequencing and by measuring virulence, growth rate, sporulation and fitness. While several lineages showed increases in virulence, we saw none of the trade-offs commonly found in obligately-killing parasites. Phenotypically similar lineages within treatments often shared multiple single-nucleotide variants, indicating strong convergent evolution. The most dramatic phenotypic changes were in timing of sporulation and spore production in vitro. We found that early sporulation led to reduced competitive fitness but could increase yield of spores on media, a trade-off characteristic of social conflict. Notably, the selection regime with strongest between-population competition and lowest genetic diversity produced the most consistent shift to early sporulation, as predicted by social evolution theory. Multi-level selection therefore revealed social interactions novel to fungi and showed that these biocontrol agents have the genomic flexibility to improve multiple traits-virulence and spore production-that are often in conflict in other parasites.

Dispersal stabilizes coupled ecological and evolutionary dynamics in a host-parasitoid system.

When ecological and evolutionary dynamics occur on comparable timescales, persistence of the ensuing eco-evolutionary dynamics requires both ecological and evolutionary stability. This unites key questions in ecology and evolution: How do species coexist, and what maintains genetic variation in a population? In this work, we investigated a host-parasitoid system in which pea aphid hosts rapidly evolve resistance to Aphidius ervi parasitoids. Field data and mathematical simulations showed that heterogeneity in parasitoid dispersal can generate variation in parasitism-mediated selection on hosts through time and space. Experiments showed how evolutionary trade-offs plus moderate host dispersal across this selection mosaic cause host-parasitoid coexistence and maintenance of genetic variation in host resistance. Our results show how dispersal can stabilize both the ecological and evolutionary components of eco-evolutionary dynamics.

High specificity of symbiont-conferred resistance in an aphid-parasitoid field community.

Host-parasite coevolution is mediated by genetic interactions between the antagonists and may lead to reciprocal adaptation. In the black bean aphid, Aphis fabae fabae, resistance to parasitoids can be conferred by the heritable bacterial endosymbiont Hamiltonella defensa. H. defensa has been shown to be variably protective against different parasitoid species, and different genotypes of the black bean aphid's main parasitoid Lysiphlebus fabarum. However, these results were obtained using haphazard combinations of laboratory-reared insect lines with different origins, making it unclear how representative they are of natural, locally (co)adapted communities. We therefore comprehensively sampled the parasitoids of a natural A. f. fabae population and measured the ability of the five most abundant species to parasitize aphids carrying the locally prevalent H. defensa haplotypes. H. defensa provided resistance only against the dominant parasitoid L. fabarum (70% of all parasitoids), but not against less abundant parasitoids, and resistance to L. fabarum acted in a genotype-specific manner (G × G interactions between H. defensa and L. fabarum). These results confirm that strong species- and genotype-specificity of symbiont-conferred resistance is indeed a hallmark of wild A. f. fabae populations, and they are consistent with symbiont-mediated adaptation of aphids to the parasitoids posing the highest risk.

Evaluating evidence for co-geography in the Anopheles-Plasmodium host-parasite system.

The often tight association between parasites and their hosts means that under certain scenarios, the evolutionary histories of the two species can become closely coupled both through time and across space. Using spatial genetic inference, we identify a potential signal of common dispersal patterns in the Anopheles gambiae and Plasmodium falciparum host-parasite system as seen through a between-species correlation of the differences between geographic sampling location and geographic location predicted from the genome. This correlation may be due to coupled dispersal dynamics between host and parasite but may also reflect statistical artifacts due to uneven spatial distribution of sampling locations. Using continuous-space population genetics simulations, we investigate the degree to which uneven distribution of sampling locations leads to bias in prediction of spatial location from genetic data and implement methods to counter this effect. We demonstrate that while algorithmic bias presents a problem in inference from spatio-genetic data, the correlation structure between A. gambiae and P. falciparum predictions cannot be attributed to spatial bias alone and is thus likely a genetic signal of co-dispersal in a host-parasite system.

The evolution of constitutively active humoral immune defenses in Drosophila populations under high parasite pressure.

Both constitutive and inducible immune mechanisms are employed by hosts for defense against infection. Constitutive immunity allows for a faster response, but it comes with an associated cost that is always present. This trade-off between speed and fitness costs leads to the theoretical prediction that constitutive immunity will be favored where parasite exposure is frequent. We selected populations of Drosophila melanogaster under high parasite pressure from the parasitoid wasp Leptopilina boulardi. With RNA sequencing, we found the evolution of resistance in these populations was associated with them developing constitutively active humoral immunity, mediated by the larval fat body. Furthermore, these evolved populations were also able to induce gene expression in response to infection to a greater level, which indicates an overall more activated humoral immune response to parasitization. The anti-parasitoid immune response also relies on the JAK/STAT signaling pathway being activated in muscles following infection, and this induced response was only seen in populations that had evolved under high parasite pressure. We found that the cytokine Upd3, which induces this JAK/STAT response, is being expressed by immature lamellocytes. Furthermore, these immune cells became constitutively present when populations evolved resistance, potentially explaining why they gained the ability to activate JAK/STAT signaling. Thus, under intense parasitism, populations evolved resistance by increasing both constitutive and induced immune defenses, and there is likely an interplay between these two forms of immunity.

Separation of evolutionary timescales in coevolving species.

Many coevolutionary processes, including host-parasite and host-symbiont interactions, involve one species or trait which evolves much faster than the other. Whether or not a coevolutionary trajectory converges depends on the relative rates of evolutionary change in the two species, and so current adaptive dynamics approaches generally either determine convergence stability by considering arbitrary (often comparable) rates of evolutionary change or else rely on necessary or sufficient conditions for convergence stability. We propose a method for determining convergence stability in the case where one species is expected to evolve much faster than the other. This requires a second separation of timescales, which assumes that the faster evolving species will reach its evolutionary equilibrium (if one exists) before a new mutation arises in the more slowly evolving species. This method, which is likely to be a reasonable approximation for many coevolving species, both provides straightforward conditions for convergence stability and is less computationally expensive than traditional analysis of coevolution models, as it reduces the trait space from a two-dimensional plane to a one-dimensional manifold. In this paper, we present the theory underlying this new separation of timescales and provide examples of how it could be used to determine coevolutionary outcomes from models.

Rapid bacteria-phage coevolution drives the emergence of multiscale networks.

Interactions between species catalyze the evolution of multiscale ecological networks, including both nested and modular elements that regulate the function of diverse communities. One common assumption is that such complex pattern formation requires spatial isolation or long evolutionary timescales. We show that multiscale network structure can evolve rapidly under simple ecological conditions without spatial structure. In just 21 days of laboratory coevolution, Escherichia coli and bacteriophage Φ21 coevolve and diversify to form elaborate cross-infection networks. By measuring ~10,000 phage-bacteria infections and testing the genetic basis of interactions, we identify the mechanisms that create each component of the multiscale pattern. Our results demonstrate how multiscale networks evolve in parasite-host systems, illustrating Darwin's idea that simple adaptive processes can generate entangled banks of ecological interactions.

Lack of host specialization despite selective host use in brood parasitic cuckoo catfish.

Host-parasite dynamics involve coevolutionary arms races, which may lead to host specialization and ensuing diversification. Our general understanding of the evolution of host specialization in brood parasites is compromised by a restricted focus on bird and insect lineages. The cuckoo catfish (Synodontis multipunctatus) is an obligate parasite of parental care of mouthbrooding cichlids in Lake Tanganyika. Given the ecological and taxonomic diversity of mouthbrooding cichlids in the lake, we hypothesized the existence of sympatric host-specific lineages in the cuckoo catfish. In a sample of 779 broods from 20 cichlid species, we found four species parasitized by cuckoo catfish (with prevalence of parasitism of 2%-18%). All parasitized cichlids were from the tribe Tropheini, maternal mouthbrooders that spawn over a substrate (rather than in open water). Phylogenetic analysis based on genomic (ddRAD sequencing) and mitochondrial (Dloop) data from cuckoo catfish embryos showed an absence of host-specific lineages. This was corroborated by analyses of genetic structure and co-ancestry matrix. Within host species, parasitism was not associated with any individual characteristic we recorded (parent size, water depth), but was costly as parasitized parents carried smaller clutches of their own offspring. We conclude that the cuckoo catfish is an intermediate generalist and discuss costs, benefits and constraints of host specialization in this species and brood parasites in general.

Leishmania genetic exchange is mediated by IgM natural antibodies.

Host factors that mediate Leishmania genetic exchange are not well defined. Here we demonstrate that natural IgM (IgMn)1-4 antibodies mediate parasite genetic exchange by inducing the transient formation of a spherical parasite clump that promotes parasite fusion and hybrid formation. We establish that IgMn from Leishmania-free animals binds to the surface of Leishmania parasites to induce significant changes in the expression of parasite transcripts and proteins. Leishmania binding to IgMn is partially lost after glycosidase treatment, although parasite surface phosphoglycans, including lipophosphoglycan, are not required for IgMn-induced parasite clumping. Notably, the transient formation of parasite clumps is essential for Leishmania hybridization in vitro. In vivo, we observed a 12-fold increase in hybrid formation in sand flies provided a second blood meal containing IgMn compared with controls. Furthermore, the generation of recombinant progeny from mating hybrids and parental lines were only observed in sand flies provided with IgMn. Both in vitro and in vivo IgM-induced Leishmania crosses resulted in full genome hybrids that show equal patterns of biparental contribution. Leishmania co-option of a host natural antibody to facilitate mating in the insect vector establishes a new paradigm of parasite-host-vector interdependence that contributes to parasite diversity and fitness by promoting genetic exchange.

Massive horizontal gene transfer and the evolution of nematomorph-driven behavioral manipulation of mantids.

To complete their life cycle, a wide range of parasites must manipulate the behavior of their hosts.1 This manipulation is a well-known example of the "extended phenotype,2" where genes in one organism have phenotypic effects on another organism. Recent studies have explored the parasite genes responsible for such manipulation of host behavior, including the potential molecular mechanisms.3,4 However, little is known about how parasites have acquired the genes involved in manipulating phylogenetically distinct hosts.4 In a fascinating example of the extended phenotype, nematomorph parasites have evolved the ability to induce their terrestrial insect hosts to enter bodies of water, where the parasite then reproduces. Here, we comprehensively analyzed nematomorphs and their mantid hosts, focusing on the transcriptomic changes associated with host manipulations and sequence similarity between host and parasite genes to test molecular mimicry. The nematomorph's transcriptome changed during host manipulation, whereas no distinct changes were found in mantids. We then discovered numerous possible host-derived genes in nematomorphs, and these genes were frequently up-regulated during host manipulation. Our findings suggest a possible general role of horizontal gene transfer (HGT) in the molecular mechanisms of host manipulation, as well as in the genome evolution of manipulative parasites. The evidence of HGT between multicellular eukaryotes remains scarce but is increasing and, therefore, elucidating its mechanisms will advance our understanding of the enduring influence of HGT on the evolution of the web of life.

Long live the host! Proteomic analysis reveals possible strategies for parasitic manipulation of its social host.

Parasites with complex life cycles often manipulate the phenotype of their intermediate hosts to increase the probability of transmission to their definitive hosts. Infection with Anomotaenia brevis, a cestode that uses Temnothorax nylanderi ants as intermediate hosts, leads to a multiple-fold extension of host lifespan and to changes in behaviour, morphology and colouration. The mechanisms behind these changes are unknown, as is whether the increased longevity is achieved through parasite manipulation. Here, we demonstrate that the parasite releases proteins into its host with functions that might explain the observed changes. These parasitic proteins make up a substantial portion of the proteome of the hosts' haemolymph, and thioredoxin peroxidase and superoxide dismutase, two antioxidants, exhibited the highest abundances among them. The largest part of the secreted proteins could not be annotated, indicating they are either novel or severely altered during recent coevolution to function in host manipulation. We also detected shifts in the hosts' proteome with infection, in particular an overabundance of vitellogenin-like A in infected ants, a protein that regulates division of labour in Temnothorax ants, which could explain the observed behavioural changes. Our results thus suggest two different strategies that might be employed by this parasite to manipulate its host: secreting proteins with immediate influence on the host's phenotype and altering the host's translational activity. Our findings highlight the intricate molecular interplay required to influence the phenotype of a host and point to potential signalling pathways and genes involved in parasite-host communication.

Cuscuta campestris fine-tunes gene expression during haustoriogenesis as an adaptation to different hosts.

The Cuscuta genus comprises obligate parasitic plants that have an unusually wide host range. Whether Cuscuta uses different infection strategies for different hosts or whether the infection strategy is mechanistically and enzymatically conserved remains unknown. To address this, we investigated molecular events during the interaction between field dodder (Cuscuta campestris) and two host species of the Solanum genus that are known to react differently to parasitic infection. We found that host gene induction, particularly of cell wall fortifying genes, coincided with a differential induction of genes for cell wall degradation in the parasite in the cultivated tomato (Solanum lycopersicum) but not in a wild relative (Solanum pennellii). This indicates that the parasite can adjust its gene expression in response to its host. This idea was supported by the increased expression of C. campestris genes encoding an endo-ß-1,4-mannanase in response to exposure of the parasite to purified mono- and polysaccharides in a host-independent infection system. Our results suggest multiple key roles of the host cell wall in determining the outcome of an infection attempt.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions.

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Different trematode parasites in the same snail host: Species-specific or shared microbiota?

The concept that microbes associated with macroorganisms evolve as a unit has swept evolutionary ecology. However, this idea is controversial due to factors such as imperfect vertical transmission of microbial lineages and high microbiome variability among conspecific individuals of the same population. Here, we tested several predictions regarding the microbiota of four trematodes (Galactosomum otepotiense, Philophthalmus attenuatus, Acanthoparyphium sp. and Maritrema novaezealandense) that parasitize the same snail host population. We predicted that each parasite species would harbour a distinct microbiota, with microbial composition similarity decreasing with increasing phylogenetic distance among parasite species. We also predicted that trematode species co-infecting the same individual host would influence each other's microbiota. We detected significant differences in alpha and beta diversity, as well as differential abundance, in the microbiota of the four trematode species. We found no evidence that phylogenetically closely related trematodes had more similar microbiota. We also uncovered indicator bacterial taxa that were significantly associated with each trematode species. Trematode species sharing the same snail host showed evidence of mostly one-sided bacterial exchanges, with the microbial community of one species approaching that of the other. We hypothesize that natural selection acting on specific microbial lineages may be important to maintain differences in horizontally acquired microbes, with vertical transmission also playing a role. In particular, one trematode species had a more consistent and diverse bacteriota than the others, potentially a result of stronger stabilizing pressures. We conclude that species-specific processes shape microbial community assembly in different trematodes exploiting the same host population.

Parasite manipulation of host phenotypes inferred from transcriptional analyses in a trematode-amphipod system.

Manipulation of host phenotypes by parasites is hypothesized to be an adaptive strategy enhancing parasite transmission across hosts and generations. Characterizing the molecular mechanisms of manipulation is important to advance our understanding of host-parasite coevolution. The trematode (Levinseniella byrdi) is known to alter the colour and behaviour of its amphipod host (Orchestia grillus) presumably increasing predation of amphipods which enhances trematode transmission through its life cycle. We sampled 24 infected and 24 uninfected amphipods from a salt marsh in Massachusetts to perform differential gene expression analysis. In addition, we constructed novel genomic tools for O. grillus including a de novo genome and transcriptome. We discovered that trematode infection results in upregulation of amphipod transcripts associated with pigmentation and detection of external stimuli, and downregulation of multiple amphipod transcripts implicated in invertebrate immune responses, such as vacuolar ATPase genes. We hypothesize that suppression of immune genes and the altered expression of genes associated with coloration and behaviour may allow the trematode to persist in the amphipod and engage in further biochemical manipulation that promotes transmission. The genomic tools and transcriptomic analyses reported provide new opportunities to discover how parasites alter diverse pathways underlying host phenotypic changes in natural populations.

Brain and antennal transcriptomes of host ants reveal potential links between behaviour and the functioning of socially parasitic colonies.

Insect social parasites are characterized by exploiting the hosts' social behaviour. Why exactly hosts direct their caring behaviour towards these parasites and their offspring remains largely unstudied. One hypothesis is that hosts do not perceive their social environment as altered and accept the parasitic colony as their own. We used the ant Leptothorax acervorum, host of the dulotic, obligate social parasite Harpagoxenus sublaevis, to shed light on molecular mechanisms underlying behavioural exploitation by contrasting tissue-specific transcriptomes in young host workers. Host pupae were experimentally (re-)introduced into fragments of their original, another conspecific, heterospecific or parasitic colony. Brain and antennal mRNA was extracted and sequenced from adult ants after they had lived in the experimental colony for at least 50 days after eclosion. The resulting transcriptomes of L. acervorum revealed that ants were indeed affected by their social environment. Host brain transcriptomes were altered by the presence of social parasites, suggesting that the parasitic environment influences brain activity, which may be linked to behavioural changes. Transcriptional activity in the antennae changed most with the presence of unrelated individuals, regardless of whether they were conspecifics or parasites. This suggests early priming of odour perception, which was further supported by sensory perception of odour as an enriched function of differentially expressed genes. Furthermore, gene expression in the antennae, but not in the brain corresponded to ant worker behaviour before sampling. Our study demonstrated that the exploitation of social behaviours by brood parasites correlates with transcriptomic alterations in the central and peripheral nervous systems.

Convergent Evolution in a Murine Intestinal Parasite Rapidly Created the TGM Family of Molecular Mimics to Suppress the Host Immune Response.

The Transforming Growth Factor-ß mimic (TGM) multigene family was recently discovered in the murine intestinal parasite Heligmosomoides polygyrus. This family was shaped by an atypical set of organismal and molecular evolutionary mechanisms along its path through the adaptive landscape. The relevant mechanisms are mimicry, convergence, exon modularity, new gene origination, and gene family neofunctionalization. We begin this review with a description of the TGM family and then address two evolutionary questions: "Why were TGM proteins needed for parasite survival" and "when did the TGM family originate"? For the former, we provide a likely answer, and for the latter, we identify multiple TGM building blocks in the ruminant intestinal parasite Haemonchus contortus. We close by identifying avenues for future investigation: new biochemical data to assign functions to more family members as well as new sequenced genomes in the Trichostrongyloidea superfamily and the Heligmosomoides genus to clarify TGM origins and expansion. Continued study of TGM proteins will generate increased knowledge of Transforming Growth Factor-ß signaling, host-parasite interactions, and metazoan evolutionary mechanisms.

Genome-Wide Analysis of Haemonchus contortus Proteases and Protease Inhibitors Using Advanced Informatics Provides Insights into Parasite Biology and Host-Parasite Interactions.

Biodiversity within the animal kingdom is associated with extensive molecular diversity. The expansion of genomic, transcriptomic and proteomic data sets for invertebrate groups and species with unique biological traits necessitates reliable in silico tools for the accurate identification and annotation of molecules and molecular groups. However, conventional tools are inadequate for lesser-known organismal groups, such as eukaryotic pathogens (parasites), so that improved approaches are urgently needed. Here, we established a combined sequence- and structure-based workflow system to harness well-curated publicly available data sets and resources to identify, classify and annotate proteases and protease inhibitors of a highly pathogenic parasitic roundworm (nematode) of global relevance, called Haemonchus contortus (barber's pole worm). This workflow performed markedly better than conventional, sequence-based classification and annotation alone and allowed the first genome-wide characterisation of protease and protease inhibitor genes and gene products in this worm. In total, we identified 790 genes encoding 860 proteases and protease inhibitors representing 83 gene families. The proteins inferred included 280 metallo-, 145 cysteine, 142 serine, 121 aspartic and 81 "mixed" proteases as well as 91 protease inhibitors, all of which had marked physicochemical diversity and inferred involvements in >400 biological processes or pathways. A detailed investigation revealed a remarkable expansion of some protease or inhibitor gene families, which are likely linked to parasitism (e.g., host-parasite interactions, immunomodulation and blood-feeding) and exhibit stage- or sex-specific transcription profiles. This investigation provides a solid foundation for detailed explorations of the structures and functions of proteases and protease inhibitors of H. contortus and related nematodes, and it could assist in the discovery of new drug or vaccine targets against infections or diseases.

Heritability and preadult survivorship costs of ectoparasite resistance in the naturally occurring Drosophila-Gamasodes mite system.

Our understanding of the evolutionary significance of ectoparasites in natural communities is limited by a paucity of information concerning the mechanisms and heritability of resistance to this ubiquitous group of organisms. Here, we report the results of artificial selection for increasing ectoparasite resistance in replicate lines of Drosophila melanogaster derived from a field-fresh population. Resistance, as ability to avoid infestation by naturally co-occurring Gamasodes queenslandicus mites, increased significantly in response to selection and realized heritability (SE) was estimated to be 0.11 (0.0090). Deployment of energetically expensive bursts of flight from the substrate was a main mechanism of host resistance that responded to selection, aligning with previously documented metabolic costs of fly behavioral defenses. Host body size, which affects parasitism rate in some fly-mite systems, was not shifted by selection. In contrast, resistant lines expressed significant reductions in larva-to-adult survivorship with increasing toxic (ammonia) stress, identifying an environmentally modulated preadult cost of resistance. Flies selected for resistance to G. queenslandicus were also more resistant to a different mite, Macrocheles subbadius, suggesting that we documented genetic variation and a pleiotropic cost of broad-spectrum behavioral immunity against ectoparasites. The results demonstrate significant evolutionary potential of resistance to an ecologically important class of parasites.